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In signal processing, tlie function of a filter is to remove unwanted parts of tlie signal, such as 
random noise, or to extract useful parts of the signal, such as the components lying within a 
certain frequency range. 
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Figure 1 : A blocl< diagram of a basic filter 



There are two main kinds of filter, analog ar 
makeup and in how they worl^. 



, They are quite different in their physical 



An analog filter uses analog electronic circuits made up from components such as resistors, 
capacitors and op amps to produce the required filtering effect. Such filter circuits are widely 
used in such applications as noise redu(^on, video signal enhancement, graphic equalisers in 
hi-fi systems, and many other areas. 

There are well-established standard techniques for designing an analog filter circuit for a given 
requirement. At all stages, the signal being filtered is an electrical voltage or current which is the 
direct analog of the physical quantity (for example, a sound or video signal or transducer output) 
involved. 

A digital filter uses a digital processor to perform numerical calculations on sampled values of 
the signal. The processor may be a general-purpose computer such as a PC, or a specialized 
DSP (Digital Signal Processor) chip. 

The analog input signal must first be sampled and digitized using an ADC (analog-to-digitai 
converter). The resulting binary numbers, representing successive sampled values of the input 
signal, are transferred to the processor, which candies out numerical calculations on them. These 
calculations typically involve multiplying the input values by constants and adding the products 
together. If necessary, the results of these calculations, which now represent sampled values of 
the filtered signal, are output through a DAC (digitai-to-anaiog converter) to convert the signal 
back to analog form. 

Note that in a digital filter, the signal is represented by a sequence of numbers, rather than a 
voltage or current 
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Figure 2: A block diagram of a basie digital filter 



Advantages of Digital Filters 

The following list gives some of ttie main advantages of digital over analog filters: 

1 . A digital filter is programmable, in other words, its operation is determined by a program 
stored in the processor's memory. This means the digital filter can easily be changed 
without affecting the circuitry (hardware). An analog filter can only be changed by 
redesigning the filter circuit. 

2. Digital filters are easily designed, tested and implemented on a general-purpose 
computer or workstation. 

3. The characteristics of analog filter circuits (particularly those containing active 
components) are subject to drift and are dependent on temperature. Digital filters do not 
suffer from these problems, and so arnt^Klfimely stable with respect both to time and 
temperature. 

4. Unlike their analog counterparts, digital filters can handle low frequency signals 
accurately. As the speed of DSP technology continues to increase, digital filters are being 
applied to high frequency signals in the RF (radio frequency) domain which, in the past, 
was the exclusive preserve of analog technology. 

5. Digital filters are very much more versatile in their ability to process signals in a variety of 
ways. This versatility includes the ability of some types of digital filter to adapt to changes 
in the charteterlste of tie signal. 

Fast DSP processors can handle complex (X)mbinations of filters in parallel or cascade (series), 
making the hardware requiremente relatively ^n^ti and compact in comparison with the 
equivalent analog circuitry. 



Operation of Digital Filters 

In the next few sections, we will develop the basic theory of the operation of digital filters. This is 
essential to an understanding of how digital filters are designed and used. First of all, we need to 
introduce a basic notation. 

Suppose the "raw" signal that is to be digitally filtered is in the form of a voltage waveform 

described by the function 

V = X (t) 

where t is time. 



fecMraaunil5f/home/20365?print 



i$mmm 



Digital Filters: An Introduction 



Page 3 of 12 



This signal is sampled at time intervals h (the sampling interval). The sampled value at time t = 
ih Is 

X, = X (ih) 

Thus the digital values transferred from the ADC to the processor can be represented by the 
s^uence 

Xqi x^, X2, x^, ... 

corresponding to the values of the signal waveform at tim^ t = 0, h, 2h, 3h, ... (where t = is the 
instant at which sampling begins). 

At time t = nh (where n is some positive integer), the values avaiiable to the processor, stored in 
memory, are 

Xq, X^, Xg, X3, ... , Xp 

Note that the sampled values x^,^,, x^^j- 3"*^ so on are not available as they haven't happened 
yet! 

The digital output from the processor to the OAQ oonsists of the sequence of values 

Yo. y1.y2.y3 - yn 

In general, the value of is calculated from the values Xq, x^, X2, X3, ... , x^. The way in which 
the y values are calculated from the x values determines the filtering action of the digital filter. 

Examples of Simple Digital Filters 

The following examples illustrate the essenfi^^glures of digital filters. 

1 . UNITY GAIN FILTER: = 

Each output value y„ is exactly the same as the corresponding input value x^: 

yo = xo 

y2 = X2 
... etc 

This is a trivial case in which the filter has no effect on the signal. 

2. SIMPLE GAIN FILTER: y^ = Kx^ (K = constant) 

This simply applies a gain factor K to each input value: 

Yo = KXq 
y-2. = KXj 

... etc 

K > 1 maizes the filter an amplifier, while < K < 1 makes it an attenuator. K < 
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corresponds to an inverting amplifier. Example 1 above is tlie special case where K = 1 . 

3. PURE DELAY FILTER: y„ = x^^ 

The output value at time t = nh is simply ttre input at time t = (n-1)h, in other words, the 
signal is delayed by time h: 

y2 = xi 

... etc 

Note that as sampling is assumed to commence at t = 0, the input value x.^ at t = -h is 
undefined . It is usual to take this (and any other values of x prior to t = 0) as zejro. 

4. TWO-TERM DIFFERENCE FILTER: y„ = \- x^.^ 

The output value at t = nh is equal to the difference betw^n the current input x^ and the 
previous input Xp.^: 

Vo = ^0 - ^-1 
yg = Xg - x^ 

yg = X3 - Xj 

... etc 

in other words, the output is the change in the input over the most recent sampling 
interval h. The effect of this filter is similar to that of an analog differentiator circuit. 

5. TWO-TERM AVERAGE FILTER: y^ = (x„ + Vi ) ' 2 

The output is the average (arithmetic rmm} ^ the current and previous input: 

yo = (Xo + Xi)/2 

y^=(Xi+Xo)/2 
y2 = (x2 + x^)/2 
yg = (X3 + X2) / 2 
... etc 

This is a simple type of low-pass filter as it tends to smooth out high-frequency variations 
in a signal. (We will look at more effeciv«s low-pass filter designs later). 

6. THREE-TERM AVERAGE FILTER: y^ = (x„ + x„.^ + x^_^) I 3 

This is similar to the previous example, with the average being taken of the current and 
two previous inputs: 

yo = (Xo + Xi + x.2)/3 
yi =(Xi +Xo + x^)/3 
y2 = (X2 + x^ + Xq) / 3 

y3 = (X3"^X2 + ^l)/3 

... etc 

As before, x.^ and x.2 are taken to be zero. 
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7. CENTRAL DIFFERENCE FILTER: = (x„ - x„.2) / 2 

This is similar in its effect to Example 4. The output is equal to half the change in the 
input signal over the current value and value two time intervals prior: 



yo 




-x.2)/2 




= 0^1 


-Xi)/2 




= {X2 


-Xo)/2 


y3 


= (X3 


-x,)/2 




etc 





Order of a Digital Filter 

The order of a digital filter can be defined as the number of previous inputs (slor^ in the 
processor's memory) used to calculate the current output. 

This is illustrBted by ttie fflters given as examples in the previous section. 

Example 1: yn = x„ 

This is a zero-order filter, since the current ou^ut depends only on the cuaent input x^ 
and not on any previous inputs. 

Example 2: y^ = Kx„ 

The order of this filter is again zero, since no previous outputs are required to give the 
current output value. 

Example 3: y„ = x„.^ 

This is a firet-order filter, as one previous input (x^.^) is required to calculate y,,. (Note that 

this filter is classed as first-order because it uses one previous input, even though the 
current input is not used). 

Example 4: y„ = x„-x„_^ 

This is again a first-order filter, since one previous input value is required to give the 
cun-ent output. 

Example 5: y„ = (x„ + x„_^) 1 2 

The order of this filter is again equal to 1 since it uses just one previous input value. 

Example 6: y„ = {x„ + x^., + x^.j) / 3 

To compute the current output y^, two previous inputs (x^,.^ and x^_2) ^i"® needed; this is 
therefore a second-order filter. 

Example 7: y„ = (x„ - x^.j) / 2 

The filter order is again 2, since the processor must store two previous inputs in order to 
compute the current output. This is unaffected by the absence of an explicit x„.^ term in 

the filter expression. 

The order of a digital filter may be any positive integer. A zero-order filter (such as those in 
Examples 1 and 2 above) is possible, but somewhat trivial, since it does not really filter the input 
signal in the accepted sense. 
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Digital Filter Coefficients 

All of the digital filter examples given in the previous section can be written in the following 
genenri forms: 



Zero ard&r. )/„ ^ 

First order: V„ = Bf^^ + a^x„., 

Second order: Yn = ao'Si * Vi * ^2^2 



Similar expressions can be developed for filters of any order. 

The constants a^, a^,^^, ... appearing in these expressions are called the filter coefficients. The 
values of these coefficients determine the characteristics of a particular filter. 

The table below gives the values of the coefficients of each of the fitters given as example in 

the previous section. 
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Recursive and Non-Recursive Filters 

For all the examples of digital filters discussed so far, the current output (y^^) is calculated solely 
from the current and previous input values (x^^, x^ x^ ^- •••)• This type of filter is said to be non- 
recursive. 

A recursive filter is one which in addition to input values also uses previous output values. 
These, like the previous input values, are stored in the processor's memory. 

The word recursive literally means "running back", and refers to 
the fact that previously-calculated output values go back into the 

calculation of the latest output. The expression for a recursive filter 
therefore contains not only terms involving the input values (x^,, x^ 

x„.2, ...) but also temns in y^.^, y^.g, ■•• 
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From this explanation, It might seem as though recursive filters 
require more calculations to be perfomned, since there are 
previous output terms in the filter expression as well as input 
terms. In fact, the reverse is usually the case. To achieve a given 
frequency response characteristic using a recursive filter generally 
requires a much lower order filter, and therefore fewer terms to be 
evaluated by the processor, than the equivalent non-recursive 
filter. This will be d^onstrated later. 



Example of a Recursive Fitter 

A simple example of a recursive digital filter is given by 

In other words, this filter determines the current output (y^) by 
adding the current input (x^) to the previous output (yn.^). 

Thus: 



and a recursive filter as an MR (or 
Infinite Impulse Response) filter. 
These terms refer to the differing 
"impulse responses* of the iwo 
types of filter 

The impulse response of a digital 
filter is the output sequence from 
the filter when a unit impulse is 
applied at its input. (A unit impulse 
is a very simple input sequence 
consisting of a single value of 1 at 
time t = 0, followed by zeros at all 
subsequent sampling instants). An 
FIR filter is one whose impulse 
response is of finite duration. An 
IIR filter is one whose impulse 
response (theoretically) continues 
forever, because the recursive 
(previous output) terms feed back 
energy into the filter input and keep 
it going. The term IIR is not very 
accurate, because the actual 
impulse responses of nearly all IIR 
filters reduce virtually to zero In a 
finite time. Nevertheless, these two 
terms are widely used. 



yi=xi+yo 
y2 = ^2 + yi 
ya = ^3 + yz 

... etc 

Note that y.., (like x,.,) is undefined, and is usually taken to be zero. 

Let us consider the effect of this filter in more detail. If in each of the above expressions we 
substitute for y^., the value given by the prmilliii@ i^ession, we get the following: 

yo = ^0 y-1 = 
yi=Xi+yo = Xi+Xo 

y2 = ^2 + yi = + Xq 

ya = ^3 y2 = ^3 + X2 + X^ + Xq 

... etc 

Thus we can see that y,^, the output at t = nh, is equal to the sum of the current input and all 

the previous inputs. This filter therefore sums or integrates the input values, and so has a similar 
effect to an analog integrator circuit. 

This example demonstrates an important and useful feature of recursive filters: the economy 
with which the output values are calculated, as compared with the equivalent non-recursive filter. 
In this example, each output is determined simply by adding twra numbers together. 

For instance, to calculate the output at time t = 10h, the recursive filter uses the expression 

yio = Xio + y9 

To achieve the same effect with a non-recursive filter (in other words, without using previous 
output values stored in memory) would entail using the expression 
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VlO = 5^10 ^9 ^8 + ^7 ^6 ^5 + + X3 + X2 + + Xp 



This would necessitate many more addition operations, as well as tlie storage of many more 
values in memory. 



Order of a Recursive (IIR) Digital Filter 

The order of a digital filter was defined earlier as the number of previous inputs which have to be 
stored in order to generate a given output. This definition is appropriate for non-recursive (FIR) 
filters, which use only the current and previous inputs to compute the current output. In the case 
of recursive filters, the definition can be extended as follows: 

The order of a recursive filter is the largest number of previous input or output 
values required to compute the current output. 

This definition can be regarded as being quite general: it applies both to FIR and IIR filters. 

For ^^mple, ttie recursive filter discussed above, given by the expression 



is classed as being of first order, because it uses one previous output value (yn.i). even though 
no previous inputs are required. 

In practice, recursive filters usually require the same number of previous inputs and outputs. 
Thus, a first-order recursive filter generally requires one previous input (x^ and one previous 

output (yn.,), while a second-order recursive filter makes use of two previous inputs (x^_^ and x„ 

2) and two previous outputs (y^.., and y^.2); and so tn, for higher orders. 

Note that a recursive (IIR) filter must, by definition, be of at least first order; a zero-order 
recursive filter is an impossibility. 

Coefficients of Recursive (IIR) Digital Filters 

From the above discussion, we can see that a recursive filter is basically like a non-recursive 
filter, with the addition of extra terms involving previous outputs (y^..,, y„.2, and so on). 

A first-order recursive filter can be written in the general form 



Note the minus sign In front of the "recursive" term b,y^.,, and the factor (1/bg) applied to all the 

coefficients. The reason for expressing the filter in this way is that it allows us to rewrite the 
expression in the following symmetrical form: 



yn = (ao>«n-^aiXn.1-biyn-l)/bo 



boyn + biy,.i = aoX„ + a^x, 



In the case of a second-order filter, the general form is 



umMm 
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An alternative "symmetrical" form of this expression is 

Note the convention that the coefficients of the inputs (the x's) are denoted by a's, while the 
coefficients of the outputs (the y's) are denoted by b's. 



The Transfer Function of a Digital Filter 

In the last section, we used two different ways of expressing the action of a digital filter: a form 
giving the output y^ directly, and a "symmetrical" form with all the output terms (y's) on one side 

and all the input terms (x's) on the other. 

in this section, we introduce what is called the transfer function of a digital filter. This is obtained 
firom the symmetrical form of the filter expression, and it allows us to describe a filter by means 

of a convenient, compact expression. The transfer function of a filter can be used to determine 
many of the characteristics of the filter, such as its frequency response. 



The Unit Delay Operator 

First of all, we must Introduce the unit delay operator, denoted by the symbol 

When applied to a sequence of digital values, this operator gives the previous value In the 
sequence. Therefore, it introduces a delay of on@ sampling interval. 

Applying the operator to an input value (say x„) giv^ the previous input (x^.^ ): 

Suppose we have an input sequence 

Xo = 5 
x,=-2 

X2 = 
X3 = 7 

x^ = 10 

Then 







= Xo = 5 




X2 


= x,=-2 




^3 


= X2 = 



and so on. Note that z'^ x^ would be x.^ which is unknown (and usually taken to be zero, as we 
have already seen). 

Similarly, applying the z'^ operator to an output gives the previous output: 
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Z"Vn=yn-1 

Applying the delay operator z"'' twice produces a delay of two sampling intervals: 

z-i(z-i x„) = rix„.i=V2 
We adopt the (fairly logical) convention 

z-iri=z-2 

in other words, the operator z"^ represents a delay of two sannpling intervals: 

This notation can be extended to delays of ttiree or more sampling intervals, the appropriate 
power of z"** being used. 

Let us now use this notation in the description of a recursive digital filter. Consider, for example, 
a general second-order filter, given in ite symmetrical form by the expression 

boVn + biy„.i + b2y„.2 = a^x^ + a^x^^ + a2X,.2 

We will make use of the following identtties: 

Substituting th^e expressions into the digital Mr giv^ 

(bo + b^z-i + b2Z-2) y„ = (ao + a^r^ + a2Z-2) x„ 

Rearranging this to give a direct relationship between the output and input for the filter, we get 

y„ / x„ = (ao + a^z-1 + a^z-^) I (b^ + b^z'^ + b2Z-2) 

This is the general form of the transfer function fer a second-order recursive (IIR) filter. 

For a first-order filter, the terms in z'^ are omitted. For filters of order higher than 2, further terms 
involving higher powers of z""" are added to both the numerator and denominator of the transfer 
function. 

A non-recursive (FIR) filter has a simpler transfer function which does not contain any 
denominator terms. The coefficient bg is regarded as being equal to 1 , and all the other b 

coefficients are zero. The transfer function of a second-order FIR filter can therefore be 
expressed in the general form 
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y„/x„ = ao + aiZ-i+a2Z-2 

Transfer Funeion Examples 

1 . The three-term average filter, defined by the expression 

Vn^^f^i + Vi +V2) 
can be written using the z'^ operator notation as 

yn = %('^ + 2-\+r2x„) 
= 1/3(1 +Z-1+Z-2)X„ 

The transfer function for the fitter Is therefore 

yn/'s, = %o+z-i+z-2) 

2. The general form of the transfer function for a first-order recursive filter can be written 

yn/Xn = (ao + aiZ-i)/(bo + biri) 
Consider, for example, the simple first-order recursive filter 

yn = Xn + yn^1 

which we discussed earlier. To derive the transfer function for this filter, we rewrite the 
filter e}^resslon using the Z'^ operator: 

0-Z-^)yn = '<n 

Reanranging gives the filter transfer fim^on as 

y„/x„ = 1/(1 -z-1) 

3. ^ a further sample, consider the second-order MR filter 

Collecting output terms on the left and terms on the right to give the "symmetrical" 
form of the filter expression, we get 

yn + 2y„.i-yn.2 = X,+2x„.^+X„.2 

Expressing this in terms of the z"^ operator gives 

(1 + 2Z-'' - z-2) y^ = (1 + 2ri + r2) x„ 
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and so the transfer functton is 

y„ / x„ = (1 + 2Z-1 + z-2) / (1 + 2z-i - r2) 
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